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Abstract 

We present new NLO sets of parton distributions in real photons based on 
a scheme invariant definition of the non-perturbative input. We compare the 
theoretical predictions with LEP data and a best fit allows us to constrain the 
parameters of the distributions. The shape of the gluon distribution is poorly 
constrained and we consider the possibility to measure it in photoproduction 
experiments. Three parametrizations which aim to take into account the 
scattering of LEP data are proposed. They are compared to other NLO 
parametrizations. 


1 Introduction 


Since the early days of QCD, the photon structure function has attracted much 
interest, and the pioneering work of Witten [I] triggered a large amount of theo¬ 
retical and experimental studies |2j. Recent developments are nicely reviewed in 
ref. HUE]. The present situation is characterized by much recent data, essentially 
accumulated by LEP experiments, by the possibility to observe the photon structure 
function in photoproduction experiments at HERA jj], and by the necessity to have 
accurate predictions for the Next Linear Collider (NLC). These three reasons justify 
an upgrading of the AFG parametrization of quark and gluon distributions in the 
photon that we proposed ten years ago f|. 

The NLO AFG parametrization was characterized by a non-perturbative input, 
defined in a factorization-scheme invariant way, and by a parameter Q q fixing the 
starting point of the (Revolution of the perturbative component. With the choice 
Ql = .5 GeV 2 (a value close to the p-rnass squared) and a non-perturbative input 
determined within the framework of the Vector Dominance Model (VDM), we found 
good agreement with data. 

Data on the photon structure function essentially determine the quark content of 
the photon. On the other hand the gluon content can be constrained in photoprod¬ 
uction reactions at HERA 0 m and the AFG gluon distribution appears to be in 
agreement with recent data on jet production |H] - However the latter lacks flexibility 
and a parametrization containing adjustable parameters should allow a better fit of 
the relevant data. In particular the VDM input used in the AFG parametrization 
rests on the 7T° structure function determined from prompt photon and Drcll-Yan 
experiments and the user is not allowed to modify this input. Moreover, the 
parametrization was only valid for Nf = 4 ; the large energies reached in collider 
experiments now require that we take into account the bottom quark contribution. 

The new AFG04 parametrization of the quark and gluon distributions in the real 
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photon is valid for Nf = 5. We work at the NLO approximation and within the mass¬ 
less, flavor changing scheme. However we keep m 2 q /Q 2 corrections (q = c,b ) in the 
direct contribution in order to have smooth thresholds when calculating F7(x, Q 2 ). 
Asymptotically, when m 2 /Q 2 goes to zero, we recover the usual MS factorization 
scheme for massless partons. The non-perturbative input, always inspired by the 
VDM approximation, has a flexible parametrization : the gluon and the sea normal¬ 
ization, as well as the gluon shape can be modified. The overall normalization of the 
non-perturbative input is also left free, and the perturbative parameter Q q can be 
varied. We study the effects on of the variation of these parameters ; constraints 
are obtained from the confrontation of the theoretical predictions with LEP data. 
As expected data on F^ do not give access to the gluon content of the photon. A 
better determination of the latter should be obtained from large-p_i_ photoproduc¬ 
tion reactions that we briefly consider. A default parametrization results from these 
studies. Other parametrizations, which reflect the scattering of LEP data, are also 
proposed. 

In section 2, we discuss the necessity to introduce a scheme-independent non- 
perturbative input. The method to reach this goal is detailed in section 3 and 
appendix A. In section 4, we present a specific non-perturbative input obtained 
from the Vector Dominance Model. Section 5 is devoted to the study of medium-Q 2 
LEP data, which allows us to constrain the parametrization of the distributions. We 
propose three different distributions which take into account the scattering of LEP 
data, and compare our best fit parametrization to the GRS TU and CJK ill} NLO 
parametrizations. Finally the gluon distribution is considered in detail in section 
6. Appendix A presents a derivation of the scheme-invariant non-perturbative for¬ 
malism, and Appendix B presents the parametrizations of the parton distributions 
available in the form of a FORTRAN code. 


2 



2 Scheme Invariant Non-Perturbative Input 

In this section we recall the method we used [5] to study the link between the 
non-perturbative and the perturbative components of the photon structure function. 
Once this link is understood, a factorization scheme invariant non-perturbative com¬ 
ponent can be defined. 

Let us start with a few definitions. The evolutions of the gluon distribution 

Nf 

G' r (x,Q 2 ), of the singlet distribution S 7 (x, Q 2 ) = E (ql(x, Q 2 ) + ql(x,Q 2 )) = 

/=i 

J2fq { f + \x,Q 2 ) and of the non-singlet distributions q^ s (x,Q 2 ) = q — E 7 /iVy (Nf 
is the number of flavors) are governed by the inhomogeneous DGLAP equations j!2| 


d log Q 2 ~ k qJr P qq ® ^ + ^ qg ® ^ 

(2.1a) 

dG 1 

dlogQ2 = k 9 + P 9q ®^ + p 99 ® G1 

(2.1b) 

dn NS 

uq f _ „NS U , D ,0, „NS 

aiogQ 2 ~ f + 

(2.2) 

where a? s = (e)/ < e 2 > -1 )/N f with < e m >= E ef/N f . 
defined by 

The convolution ® is 

P0q= 

Jx Z \zj 

(2.3) 

The homogeneous (P%j) splitting functions were calculated 

in ref. |12, 33] th e 

NLO approximation. The inhomogeneous splitting functions 


/, _ a jJo) , a a s(Q 2 ) , (l ) 
kq ~2ir^ + 2vr 2vr kq 

(2.4) 

a a s (Q 2 ) (1) 

9 2 tt 2 tt 9 

(2.5) 
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may be derived from the Pj 3 and are given in refs. . T51 US] ; the expression for the 
LO splitting function is 2Nf < e 2 > 3[x 2 + (1 — x) 2 ]*. 

In terms of the parton distributions, the photon structure function is written 


Tl(x, Q 2 ) = F2{x, Q 2 )/x = E e / 4 +) ® C q + G 7 ® C a + C 7 . (2.6) 

/ 

The Wilson coefficients C q and C g may be found in ref. |2XJ, and the direct term 
C-y, in the MS scheme, is given by EHEU 


C 1 = 2 ^ 2 E e / 3 + (! - ^) 2 ) ln 


1 — x 


x 


+ 8x(l — x) — 1 


(2.7) 


The physical quantity is factorization scheme independent. This means that 
it does not depend on the procedure (the factorization scheme) used to define the 
NLO splitting function (n > 1) and k\ n \ and the function C q , C g and C T 
This is however true only if these functions were calculated to all orders in a s . If 
the truncated series (Q and m are used, the photon structure function is still 
scheme independent, but only at order 0(ot Q a ). 

Let us consider, for the sake of simplicity, the non-singlet eq. m- Its solution 
can be written, for moments of the quark distribution qf s (n) = /q 1 dxx n ~ 1 q 1 ^ s (x, Q 2 ), 
as follows : 


?f s (n) = af s j Q2 ^k q (n)e!S (2 . 8 ) 

For small values of k 2 , the perturbative approach is no longer valid. Let us assume 
that we can use this expression for k 2 > Qq ; we then define the perturbative 
(Anomalous [I]) component 


-PNs(n) 


<1 N a S n(Q\QI) = o- 


NS 




Pns 


(2.9) 


*We do not consider NNLO corrections. Therefore our parametrizations are consistent with 
the NLO calculations of large-pr photoproduction cross sections min]. Expressions of NNLO 
corrections and discussions of their importance may be found in GBJDJEOi. 
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(we have dropped the indices / and n). 

For k 2 smaller than Ql, we are in the realm of non-perturbative QCD and we 
write the corresponding hadronic contribution (which behaves like a hadron struc¬ 
ture function and is discussed in detail in appendix A) 


q N H S {Q\Ql) = q N H ^Q 


NS/ 


rQ 2 dk’‘ 
Jq 2 


e jQ o 


Pns 


( 2 . 10 ) 


the total non-singlet distribution being the sum of the Anomalous and Hadronic 
component 


q NS m = q n a s n{Q\Q l) + 4 s (Q 2 ,Ql) • (2.11) 

Actually with EH we have written the general solution of the inhomogeneous 
equation El, the only assumption being that the scale Ql allows us to define 
a perturbative and a non-perturbative component. However this way of defining 
a non-perturbative component is too naive and factorization scheme dependent. 
Indeed let us consider the contribution of the HO inhomogeneous kernel 
to the anomalous component El 


q^(Q‘ 2 :Ql) = 


' a s m 

. a a {Ql . 


-2 P$/0o 


a ^ a NS _|_ 


2vr p(P) 


where we used the running coupling constant dehned by 


( 2 . 12 ) 


da s (Q- 


W =- Q '(4^ + (S) A) • (2-13) 

When similar expressions for q^ are introduced in (12.61) . one obtains for JQ a 
contribution proportional to a® (besides the Leading Logarithm contribution) 


n ~ - 


a 




< e 4 > 


27r pj°) < e 2 > 


+ a 


( 2 . 14 ) 
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which is factorization scheme independent, and a contribution which verifies the 
homogeneous LO DGLAP equation 






h w e / 

27r p® TV/ < e > 2 



/'o,(« 2 )V 2P " /ft 

ww 


(2.15) 


which must be also scheme independent. Now it is clear from (ITT5D that q^j is not 
scheme invariant with respect to the “photon factorization” scheme which defines 
the inhomogeneous kernel k' q l> . Therefore it cannot be, for instance, the same in the 
MS scheme or the DIS 7 scheme defined in m- Of course, q £ is also non-invariant 
with respect to the usual hadronic factorization scheme which defines P/£ Thus 
the assumption that the hadronic input could be described by a VDM-type input is 
clearly too naive. It may however be true in a specific factorization scheme and we 
explore this possibility in the next section and in appendix A. 


3 The Non-Perturbative Input at Lowest Order 

In order to better understand the content of qii(Qo), let us consider the lowest 
order contribution to coming from the imaginary part of the box diagram, Fig. 1, 
which shows how the virtual photon q probes the quark content of the real photon 
p. The lower part G(fc,p)/(— k 2 ) (it includes the quark propagators) represents the 
coupling of the real photon to a qq pair and includes non-perturbative effects. 

Actually our only assumption is that G(k,p) tends to the pointlike term for large 

\k 2 l 


Jim G(k,p) = G p (k,p) ~ S ((p- kf) [z 2 + (1 - z) 2 


|fc 2 |>A : 


(3.1) 


where z is the fraction of the longitudinal p momentum carried away by k. When 
k 2 goes to zero, G(k, p)/(—k 2 ) must be integrable, because Tk is a physical finite 


quantity. This means that we must have lim G(k,p ) ~ (— k 2 ) a with a > 0. 

k 2 —>0 
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Figure 1: The box diagram 

We make the pointlike content of G(k,p ) explicit by writing 


G(k,p) G F (k,p) G(k,p)-0(\k 2 \-Ql)G F (k,p) _ e(Ql-\k 2 \)G F (k,p) 

_p _ J/2 2 

(3.2) 

The first term on the RHS of (E 2 D, without cut on k 2 , corresponds to the pertur¬ 
bative expression of the box diagram. 

Its contribution, in the collincar approximation, is easily calculated |B| and is 
equal to 


n 


3e1 


a 


7r 


(.x 2 + (1 - x) 2 ) - 


£ \/l 2 


+ ( x 2 + (1 — x) 2 ) ln(l — x) + 2 x(l — x) 


= 3e) ^ + In (x 2 + (1 - x) 2 ) + C(£ c (3.3) 

in which we add a quark / and an antiquark / contribution. We define Q 2 = —q 2 
and x = Q 2 /‘Ip-q. Note that we took the upper bound of the integral over \k 2 \ equal 
to Q 2 . (The actual bound is Q 2 /x, but this x-dependence is beyond the collincar 
approximation). 

This expression for the box diagram has been obtained with the dimensional reg¬ 
ularization (| = 7 — 7 e + In 47 t); it is the one used to define the MS factorization 
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scheme which consists in subtracting the term proportional to (Q 2 / p 2 )~ £ / e. This 
procedure defines the scheme-dependent direct term C/ c in the collincar approxi¬ 
mation (or C 1 given in (EH) when we take into account the non-collinear terms). 

being a physical quantity, it cannot contain the 1/e pole, and it is here 
that the third term of the RHS of (E3) plays its part. We obtain from this last term 

^ = -3eJ l + In ) (V + (1 - *?) - Cl . (3.4) 

This term has no anomalous In Q 2 behavior. Actually it is independent of Q 2 , 
when QCD is not switched on. When the part of JF/’ c proportional to (—| + In -j^) 
is added to Q, the 1/e poles cancel each other and In Q 2 / p 2 is changed into 


, Q 2 , Q 2 

111 m = 111 v 




4(/r 


«,(Q 2 

MQl 


(3.5) 


A 2 Poa s (Q 2 ) V 

This Q 2 -dependence corresponds to the LO part of m wit fi p$ = o. 

Let us now consider the second term of (jH. The 0-function cuts the 1/A; 2 
perturbative behavior of this contribution. The integration over k 2 is therefore 
controlled by the non-perturbative behavior of G(A;,p)/(— k 2 ) and we obtain a result 
which does not depend on Q 2 . The value of Ql must of course be chosen such 
as G NP (p,k,Ql) = G(k,p ) — 0(|A; 2 | — Ql)G p (k,p) represents the non-perturbative 
input. For instance an overly large value of Ql would conduce to a perturbative tail 
in G NP (p, k, Q%). 

We now define the non-perturbative quark content of the real photon by 


cTk 5 z 


k ■ n\ G NP (p : k, Ql) 


p ■ n 


-k 2 ) 


= q 


NP, 


Z,Ql ) 


(3.6) 


(n is a light-cone vector such as k ■ n rs./ k° + At). With (13.61) we have defined a 
non-perturbative input (if Ql is correctly chosen) which is invariant with respect to 
the photon factorization scheme. Indeed it does not depend on the regularization 


used to calculate (Ed nor on the subtraction defining the MS scheme. When the 







QCD evolution is switched on (all order QCD expressions are discussed in appendix 
A), both q NP and C^ c acquire a hadronic Q 2 -dependence and we obtain a hadronic 
contribution (which behaves like a hadronic structure function) to J-% 


qfm = 


'Oi s (Q 2 


-2P$/p 0 


V(OS) - ^jr 


(3.7) 


This hadronic contribution is scheme dependent because of the presence of C^ c , 
but qf P (Qo) is not. Therefore, in the MS factorization scheme, the hadronic input 
is given by expression (ED, and, at Q 2 = Q%, we have 




N f 

2 \X,QI) = c 7 (x) + 5Z [e 

/=i 


(<tf p (Qo) + <if P (Ql)) 


C L 


(3.8) 


In the above expression, we only studied the part of T2 associated with the 
quark contributions. (Nor did we write the convolution with the Wilson coefficient). 
Similar considerations applied to the gluon distribution lead to modifications of the 
input starting at order 0(a s ). The NNLO corrections are not considered in this 
paper. 

These results are different from those obtained by the authors of ref. nnu2smni 
who work in a factorization scheme called DIS 7 in which 


a < e 4 > 
2n < e 2 > 


k\ l \DIS y ) 


a < e 4 > 
27r < e 2 > 


kf>(MS) - C 7 P<°> , 


so that Cj(DISj) = 0. In this case, the structure function is written 


(3.9) 


N f 

•^~2 ( x i Qo) = 5Z e f (Vf^DIS-yiQo) + 9/^D/5 7 (Qo)) (3.10) 

/=1 

where q^ois^Q 2 ) non-perturbative input in this particular factorization scheme. 

Let us finish this discussion by emphasizing the fact that the parton distributions 
defined by ETD and ED are universal (independent of the particular reaction 
studied here, namely the DIS on a real photon). Of course they are factorization 
scheme dependent and here we work in the MS scheme. 
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4 The Vector Dominance Model 


The non-perturbative contribution defined in (EH) is not known. This is why we 
could proceed as in the pure hadronic case by defining a parameter-dependent input 
and by determining the parameters by a fit to data. Here we prefer to follow another 
path and to try to constrain the non-perturbative input by assuming that it can be 
described by the quark and gluon distributions in Vector Mesons. This assumption, 
the Vector Dominance Model, is known to work well in the non-perturbative domain 
and to correctly describe how photons couple to quarks. We used this assumption 
in our preceding paper ||, which led to the AFG parametrization. Here we keep 
this approach, but we make it more flexible by varying the non-perturbative nor¬ 
malization of the gluon and sea quarks. We also consider modifications of the gluon 
x-shape. 

In ref. p] we considered the photon as a coherent superposition of vector mesons 

7= T 2 ( p+ 1" t*) = 9 (r s - l dl ~ l sS ) (41) 

with a coupling constant g determined from the <T ta t{'yp) and a tot (jrp) cross sections 

g 2 ~ a . (4.2) 

Assuming that the parton distributions in the qq “bound states” of EU) are similar 
to those of the pion, observed in Drcll-Yan and direct photon reactions PJ, we can 
write 

4 

U valence ( X , Q 2 ) = 9* g U valence ( 4 ‘ 3a ) 

Ulea( x , Q 2 ) = ^ + l + U ^ a = 3 2 \ U 1ea ( 4 -3b) 

G 7 (x, Q 2 ) = g 2 ^g n {x, Q 2 ) (4.3c) 
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and so on for the parton distributions of the non-perturbative component of the real 
photon (we assume a SU(3) flavor symmetry). 

This rough approach leads to a reasonable agreement with data {3j. Here we 
would like to make it more flexible. Let us start from the parametrization of the 
pion structure function at Q 2 = 2 GeV 2 taken from ref. [H] and used in ref. [HJ 

xKalence = C v x P2 ( 1 - x) p 3 (p 2 = -48 , p 3 = .85) (4.4a) 

xq: ea = C 9 (l - xr (Ps = 7.5 , = 1.2) (4.4b) 

xg n = C g ( 1 - x) P10 (pio = 1.9) (4.4c) 

with C v = 1/B(p2,l + p 3 ) ( B(x,y ) is the beta-function) and C g = (1 +pio)(l — 
1 +ps — l+ra +po ) = -447(1 + pw), C g being determined in such a way that 

[ dxx [2 ul alence + q* ea + g*] = l. (4.5) 

Jo 

This input is in good agreement with ^(x, Q 2 ) data which mainly constrain the 
quark distributions. Therefore we leave the quark distributions fixed and we take 
Pio as a free parameter. As we shall see below, J TJ is not sensitive to variations of 
Pio, a parameter which should be constrained by photoproduction data. We also 
leave some freedom in the normalization of the distributions. First of all the overall 
normalization is allowed to vary around the value fixed in (Q and we write 

g 2 = c np ■ a . (4.6) 

Then we also consider the possibility of having a different coupling of the pho¬ 
ton to the valence distributions and to the sea and gluon distributions (this extra 
coupling could proceed through a quark loop). We parametrize this possibility by a 
modification of C s and C g 
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(4.7a) 


C s 


C„ 


1.2 


C g — C mom ■ 0.447 (1 +p 10 ) , (4.7b) 

the default value being C mom = 1.0. 

Let us end this section by discussing another input, the quark masses. Threshold 
effects due to the charm quark may be important in J^ix, Q 2 ) at large x. To 
study this problem, let us again consider the box diagram and the massive quark 
contribution to J 7 ]■ Dropping all inessential factors and neglecting terms of order 
0(m 2 /Q 2 ), one obtains iTi 

Flm ~ ~P){(x 2 + ( X - ^) 2 ) ln £7 + (x 2 + (1 — ^) 2 ) In + 

+ (8a;(l-x)- 1)/T^} , (4.8) 

with s = (p + q ) 2 and (3 = Am 2 /s. Subtracting the term proportional to log we 

•i TYl q 

define a massive direct term 


1 a. 


C^(x,m q ) = e):—30(1 
7r 


P) (x 2 + (1 


+ (8x(l - x) - l)^/l - f3 


(4.9) 


which has the massless limit El when fi — 0. 

Then the term proportional to In £2 is replaced by the one generated by the 
massless evolution equations m and El with the boundary conditon q(x, Q 2 = 
m 2 ) = 0. This evolution exactly reproduces, at the lowest order in a s , the In £ 2 - 
term of (El. Therefore, close to the threshold Q 2 = rn 2 , of expression El 
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reproduces the behavior (Oil of the box diagram contribution. Let us also note that 
we keep in C 1 {x,m q ) the terms of order 0(m 2 /Q 2 ) when Q 2 ~ m 2 . 

The effect of the massive direct term (El is important when x goes to x t h = 
1/(1 + 4 m 2 /Q 2 ). One then gets In {l+ V^N) - ~ 2a/1 — f3 and given by (14.81) goes 
to zero, whereas the use of the massless limit (12.71) (without cut on x) leads to a 
negative contribution when x goes to 1. 

ffowever one must keep in mind that in most applications, we are far from the 
threshold and the massless evolution of the charm distribution is a good approxi¬ 
mation which allows us to take into account the effects of the QCD evolution, not 
present in EHD- But this is not true for the bottom distribution as long as m 2 /Q 2 is 
large. The distributions presented in this paper are obtained by solving eq. EH) and 
(12.21) with Nf — 3 for Ql < Q 2 < m 2 c) Nf = 4 for m 2 < Q 2 < ml (m c = 1.41 GeV) 
and Nf = 5 for rri'l < Q 2 (m b = 4.5 GeV). 

5 Analysis of LEP data 

In this section we analyze data on JFJ in the light of the parametrization dis¬ 
cussed in sections 3 and 4. Once we assume that the non-perturbative input can be 
determined within the framework of the Vector Dominance Model as explained in 
the preceding section, the number of free parameters is considerably reduced. JF 2 
is barely sensitive to the gluon distribution parameters C mom and pio that we shall 
discuss in relation to photoproduction reactions and, for the time being, we keep 
these parameters equal to their default values C morn = 1.0 and p 0 = 1.9. Therefore 
only two free parameters remain, C np which fixes the overall normalization of the 
non perturbative input (expression (I4.0j) i and Ql which fixes the boundary between 
the perturbative and the non-perturbative model. 

If data were precise enough, it should be possible to constrain C np and Ql sepa¬ 
rately. In fact the VDM contribution decreases rapidly with x and the perturbative 
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contribution is dominant at large x. Therefore in this kinematical domain, it should 
be possible to determine by the means of medium Q 2 data (indeed for overly 
large Q 2 , JFj is 110 longer sensitive to Qq). At small values of x on the contrary the 
non-perturbative input is large and data should constrain C np . As we shall see this 
ideal situation is not realized, because the data existing at large values of x are very 
poor, and we are led to fit low-x data where C np and are correlated. 

First let us concentrate on LEP data at small and medium Q 2 , an overall com¬ 
parison with all existing data will be conducted at the end of this section. The data 
that we analyze are in the range 3.7 < Q 2 < 17.3 GeV 2 and belong to the four LEP 
experiments. ALEPH [2SJ LEP2 data at 17.3 GeV 2 , DELPHI j27j LEP1 data at 
12.7 GeV 2 and 5.2 GeV 2 , L3 [2B] LEP2 data at 15.3 GeV 2 , and OPAL j2H| LEP1 
( Q 2 = 3.7 GeV 2 ) and LEP2 ( Q 2 = 10.7 GeV 2 ) data. The comparison between NLO 
theory and data is done in Figs. 2 and 3 where we show the theoretical curves 
obtained for Q q = .3 GeV 2 and for Ql = 1.0 GeV 2 with C np = 1. We work in the 
MS scheme and use A^^ = 300 MeV, a value which is in agreement with the world 
average and a determination obtained by fitting photon structure function data [22]. 
The coupling a s (Q 2 ) is obtained by exactly solving eq. (12 .1 311 . 

In Fig. 2, we see that large-x data are poor. Either they have large error bars 
(ALEPH, L3), or they correspond to large x-bins (DELPHI, OPAL)L Therefore 
they cannot accurately constrain Qq and do not allow us to check the predicted x- 
dependence of F^x^Q 2 ). On the other hand, all large-x data points (x > .2) favor 
a value of Ql close to 1 GeV 2 ; the choice Q q = .3 GeV 2 leads to predictions well 
above the experimental points. This is also true for smaller values of x where the 
VDM contribution is large. However in this x-region we expect a strong correlation 
between Q q and C np . 

' The errors are the total errors (when given by the experiments), or the linear sum of the syst. 
and stat. errors (OPAL). We also investigated the effect of quadratically summed errors in the fit 
performed below. Our best fit parameters change by less than 10 %. Correlations between data 
points are not taken into account. 
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Figure 2: Theory confronted with ALEPH pSj, DELPHI j2Zj, L3 [2SJ and OPAL 
data in the range 10 GeV 2 < <Q 2 > <18 GeV 2 . 


The low-Q 2 data of Fig. 3 lead us to the same conclusion. However whereas the 
agreement between OPAL data and theory is good, suggesting a slightly smaller 
contribution of the VDM component, DELPHI data on the contrary suggest a very 
large suppression, by some 30 %, of the VDM component. It is worth noting that 
the DELPHI points at x = .05 and x = .17 are below the corresponding OPAL 
points, although the Q 2 value is larger. Other data at low Q 2 do not clarify the 
situation. PLUTO data IBSj are close to the OPAL points, but the error bars are 
very large. L3 data [3T] ( Q 2 = 5 GeV 2 , not shown) are noticeably above OPAL data 
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Figure 3: Theory confronted with low-Q 2 data from OPAL, DELPHI and PLUTO 
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at small-a;, x ~ 1CU 2 , and TPC-2y [22] data at Q 2 = 5.1 GeV 2 (not shown) have 
large error bars. 

In order to determine the values of Ql and C np , we proceed to a best fit of 
the LEP data displayed in Figs. 2 and 3. The theoretical values are calculated 
with all parameters kept fixed, except Ql and C np , and they are averaged over the 
corresponding experimental x-bins. For a given value of Qq, we look at the value 
of C np which minimizes the y 2 value, and we obtain the results shown in Fig. 4 
for the ALEPH and DELPHI experiments. At the minimum of the curves, the 
corresponding values of the non-perturbative normalization are respectively C np =.60 
and C np = 1.05. 

We obtain very different shapes of the y 2 -curves. In fact, L3 and OPAL data 
lead to y 2 -curves very similar to the one displayed in Fig. 4 for ALEPH. Only 
DELPHI data lead to a false minimum at Qq=1.5 GeV 2 and C np — 1.05; there is no 
true minimum for Q g > .2 GeV 2 . It is easy to find the reason for this behavior. 
The first point at small-x of DELPHI data is high with respect to the next point. 
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Figure 4: x 2 -values as a function of Q g 

This configuration, associated with small errors, drives the fit to small and C np 
values and reproduces the steep slope of F^(x, Q 2 ) at small-x. Such an effect is not 
present in the other sets of data. The small errors of the DELPHI data give an 
important weight to this experiment in a fit of all the data sets shown in Fig. 2 and 
3. Therefore when we perform such a fit, we find (Fig. 5 (left)) a result similar to 
the DELPHI fit of Fig. 4. 

As explained in the beginning of this section, we selected medium and low-Q 2 
LEP data in order to constrain the value of Qq. This led us to consider the LEP1 
data from DELPHI. However in the LEP2 data from the same experiment, the small- 
x behavior observed in the LEP1 data is less marked. The DELPHI collaboration, 
using different Monte Carlo generators, noticed a strong model dependence of the 
resulting data on F^x^Q' 2 ). For instance the small-x data at Q 2 = 19 GeV 2 show 
a sizeable dependence on these generators and moreover, two (over three) of the 
values obtained for F^x^Q 2 ) are smaller than the one obtained at LEP1 ( Q 2 = 
12.7 Gev 2 ). Therefore, in order to take this scattering into account, we change the 
errors in the first x-bin of the LEP1 DELPHI preliminary data by a factor two. After 
this modification, we obtain the y 2 -curve shown in Fig. 5 (right). The minimum 
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Figure 5: y 2 -values of the overall fit. Without modification of the DELPHI errors 
(left) and with modification (right) 


value of the curve at Ql ~ .7 GeV 2 is obtained for C np ~ .78, and the y 2 by degree 
of freedom is x% — 1-03. 

This result demonstrates the good agreement between our theoretical input (EED 
and LEP data. However, as already discussed, these data do not well constrain the 
large-x behaviour of and the good y 2 shown in Fig. 5 cannot be seen as a 
confirmation of our large-x theoretical expression. On the other hand we must note 
that some LEP data are marginally compatible*. This is shown in Fig. 6 in which we 
display the Ay 2 = 1 contour in the Ql and C np plane. The best fits to the individual 
data sets are also exhibited and are scattered outside the contour. For instance the 
L3 data compared to theory calculated with the overall best fit parameters (Q% = .7, 
C np = .78) lead to y 2 = 4.4 which must be compared to the value obtained at the 
L3 best fit point (.7,1.1) y 2 = .85. For OPAL (10.7) we obtain y 2 = 1.76 compared 
to x 2 = .072. 

In order to partially take into account this scattering, we provide three parametriza- 
tions compatible with the contour of Fig. 6, corresponding to three different values 

*Note also that some data are not corrected for the limit P 2 —> 0 of the target photon virtuality. 
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Figure 6: The Ay 2 = 1 contour in the Q jj — C np plane with the individual best fits 
of the LEP data. The DELPHI (12.7) point (Qjj — 1-5 GeV 2 , C np = 1.05) is outside 
the figure and DELPHI (3.7) has no minimum for Qq < 1.9 GeV 2 . 

of Qq : Ql = .34 GeV 2 , Ql = .70 GeV 2 (best-fit parametrization) and Ql = 
.97 GeV 2 . We use the notation^ AFG04 (Qo,C np ,pio) and the short-hand nota¬ 
tions AFG04_BF = AFG04(.7, .78, 1.9), AFG04_LW = AFG04(.34, .6, 1.9) and 
AFG04J4G = AFG04(.97, .84, 1.9) for the parametrizations presented in Appendix 
B. 

We would obtain similar results from the world data on F%, but with a larger 
scattering of the various experiments. This can be observed from Figs. 7, 8 and 9 
in which we compare our best fit prediction with these data H3-HE]. Whereas the 
overall agreement is good, some data are clearly outside the general trend represented 
by the best fit. Therefore a general fit leading to a single parametrization has not a 
clear meaning. It is why we prefer to “frame” the data by several parametrizations, 

as we did after the analysis of recent low and medium Q 2 data sensitive to Qq and 
n 

w 77/p* 

§The parametrization AFG04(.5, 1., 1.9) has been used in ref. [UJ under the name AFG02 and 
in ref. 7 under the name AFG04. 
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Figure 7: Comparison between the best fit structure function and the world data 
with 3.7 GeV 2 < <Q 2 > < 11 GeV 2 . 

Let us now compare, in Fig. 10, the Low-Qq and High-Q^ AFG04 parametriza- 
tions to the best-fit parametrization. Let us consider the figure at the left. At small 
values of x where the non-perturbative component is large, the ratio partly reflects 
the values of C np used in AFG04J3F and AFG04_LW. At very small values of x, the 
perturbative contribution becomes dominant and the ratio reflects the effect of the 
(Revolution which is larger for AFG04_LW. At x ~ .5, the non-perturbative contri¬ 
bution is smaller and the ratio also reflects the effect of the (Revolution of the per- 
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X 


Figure 8: Comparison between the best fit structure function and the world data 
with 11 GeV 2 < <Q 2 > < 29 GeV 2 . 

turbative component. For large-a; values, the w-quark distribution contains a term 
proportionnal to ln(l — x). Adding the contribution from k ( q l> (12.121) to the Cl f c . con¬ 
tribution (18.71) . we obtain a term proportionnal to -^(1 + (a s (Q 2 ) I a s (yQl))^ 2P ^ q) ^°) 

n 

which is large for a small evolution {P qq (ji) being negative at large n ). The u-quark 
ratio reflects this behavior. Finally we note that the variations in the distribution 
functions never exceed 10 %. 

We end this section by a few words on other NLO parametrizations. Let us 
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Figure 9: Comparison between the best fit structure function and the world data 
with 29 GeV 2 < <Q 2 > < 700 GeV 2 . 

start with the AFG one. In fact the AFG parametrization is very close to the 
parametrization AFG04(.5, 1., 1.9). The only differences come from the values of 
Ails (300 MeV in AFG04 and 200 MeV in AFG) and from the absence of bottom 
quark distribution in AFG. The comparison is done at Q 2 = 50 GeV 2 , a scale which 
is in the range of those used in large-pj_ photoproduction. The ratio of the best-fit 
distribution AFG04_BF to the AFG distribution is displayed in Fig. 11. The smaller 
normalization of the non-perturbative component of AFG04J3F compared to the 
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Figure 10: The Low-Qq (left) and the High-Qg (right) parametrizations compared 
to the best-fit parametrization. 

one used in AFG explains the pattern at small-x values, ffowever, at very small 
values of x , the inhomogeneous kernels k q 1 ' 1 and k (12.41 and EB have a singular 
behavior ~ ln 2 x and k ( g l> ~ -f) and the behavior of the perturbative gluon 

and sea components reflect a faster evolution of AFG04 with Q 2 due to the larger 
value of Ajfg. At large-a; values the perturbative contributions are important and 
the ratio of the w-quark distribution is close to the ratio r = a s (AFG) / a s (AFG0A) 
at Q 2 = 50. GeV 2 , namely r ~ .9. The difference in the predictions for F^ are 
illustrated in the figure on the right. It is not very large, and could be distinguishable 
only at large values of x where data are poor. Comparisons between AFG04_BF 
and AFG predictions are given for OPAL (3.7), OPAL (10.7) and ALEPH (67.2) in 
Figs. 7 and 9. 

A comparison with the GRS ra parametrization can also be performed on the 
basis of Figs. 7 and 8. In the second reference of mu. a comparison is made between 
OPAL data and the GRS predictions which can hardly be distinguished from our 
best fit curves at Q 2 = 3.7, 10.7 and 17.8 GeV 2 . However note that the large-a: 
behaviors of F% (outside the data range) are quite different between the GRS and the 
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AFG04J3F parametrizations. This behavior comes from the different factorization 
schemes used (MS versus DIS 1 (123), associated with different non perturbative 
inputs (compare (13.81) and (13.101) ). 

In a recent publication, the authors of ref. Em proposed new 5 flavor NLO 
parametrizations and did detailed comparisons with world data, and other (AFG 
and GRS) parametrizations. On the basis of Figs. 7, 8 and 9, and of similar figures 
in ref. mi, we can easily observe several differences between the parametrizations. 
At small x (x < 10” 2 ) and for Q 2 > 5 GeV 2 the parametrizations of ref. TT] are 
higher than AFG04_BF (which is very close to GRS). This trend increases with Q 2 . 
At medium x in the charm threshold region and at larger values of x (x > .6), the 
differences between AFG04J3F, and the FFNScjkI and CJK NLO parametrizations 
of ref. Em are also noteworthy. Here also the origin of this difference is the different 
non perturbative inputs associated with the MS scheme (expression (13.81) ) and the 
DIS 7 scheme used in ref. mi- In all cases data are not accurate enough to enable 
to distinguish between the models. 



Figure 11: The ratios AFG04J3F/AFG at Q 2 = 50. GeV 2 (left) and the structure 
functions F^(x, Q 2 ) (right). 
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6 The gluon content of the photon 


In this section we study other possible options for the parametrization of the 
parton distributions in the photon. First we study a modification of the sea quark 
distributions. Then we investigate the effect of changing the normalization of the 
sea quark and gluon distributions. And finally we modify the large-x shape of the 
gluon distribution. As we shall see, these modifications are poorly constrained by 
F 2 data, but some of them could be visible in photoproduction experiments. 

The small -x behavior of the sea quark distribution (I4.4bjl that we used until now 
is less steep than those of some recent parton distributions in the proton gang ED| . 
For instance we have x(u + d) ~ .061/ar 3 at x < 10~ 3 for the CTEQ6M distribution 
at Ql = 1.69 GeV 2 . I 11 order to explore the effect of such a steep behavior, we 
modify our sea distribution (14.4bjl . while keeping fixed the momentum carried by 
the sea quarks (Jq dx x gj ea (x, Ql = 2 GeV 2 ) = .14) 

Wsea = - X) 7 ' 5 /X- 3 . ( 6 . 1 ) 

This ansatz corresponds to quite a large sea at small-rc, larger by more than a factor 
2 than the corresponding CTEQ6M parametrization for the proton. Therefore m 
must be considered as an extreme parametrization. 

The resulting y 2 is less satisfactory than the one obtained in the preceding sec¬ 
tion. With the exception of DELPHI, the individual y 2 deteriorate. The total best 
fit now corresponds to y 2 = 34.3 (without DELPHI-errors modification), close to 
the value of section 5. If we modify the errors, we obtain y 2 = 31.4, instead of 25.8 
in section 5. From these results, we see that there is no compelling reason to modify 
the small-x behavior of the sea distributions we used in the preceding section. 

Let us now study the effects of the parameter C morn which modifies the normal¬ 
ization of the sea quark and gluon distributions according to TT~T7T1 rTTbli (at small-x 
values, sea quark and gluon distributions are strongly coupled and we modify their 
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Figure 12: Correlation between the parameters C mom and C np . The numbers are 
the x 2 -values. 

normalizations by the same parameter C mom ). If we keep and C np at the best 
£t values established in section 5, C mom is quite constrained by data and the 
X 2 - value varies by less than one unit if we stay in the domain .92 < C mom < 1.08. 
However it is clear that we can partially compensate the C mom variations by also 
varying C np . Keeping Ql = .7 GeV 2 , we hrst observe a strong correlation between 
C np and C mom (Fig. 12). By playing with the values of C mom and C np , for instance, 
we can enhance the importance of the valence compared to the sea qnark. But it 
is unlikely that photoproduction experiment could better constrain C np and C mom 
and we do not pursue this study in detail. 




Figure 13: y 2 variation as a function of the gluon parameter pu). 
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Finally we consider the modification of the gluon distribution and we vary the 
parameter pio of expression m . As expected is not sensitive to the gluon 
distribution and LEP data do not constrain the value of pi 0 . We display in Fig. 13 
the dependence of the y 2 on p 10 which is weak (Q(j and C np being fixed at the best 
fit values). In a large range in pio, y 2 varies by less than one unit. In Fig. 14 we 
show the behaviour of the distributions obtained with pio = 1.0 (hard gluon) and 
P 10 = 4.0 (soft gluon) at Q 2 = 50. GeV 2 , the other parameters being kept fixed at 
the best-fit values. The behavior at small values of x is due to the normalization 
factor C g (THcl) ; at large values of x the non-perturbative inputs vanish and the 
ratios go to one. 



Figure 14: Ratios of the gluon distributions with p 10 = 4.0 (soft gluon) and p 10 = 1.0 
(hard gluon) to the best-fit gluon at Q 2 = 50. GeV 2 . 

On the other hand photoproduction reactions [El are sensitive to the gluon 
distribution since an initial gluon can interact with a parton from the initial proton 
producing two large-p^ jets in the final state. Jet production in photon-photon 
collisions is also a reaction allowing to observe the gluon distribution in the range 
0 < x < .4 J52J. A particularly interesting reaction is the photoproduction of large- 
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p l photon and jet, which has been studied in ref. [TJ in detail. The interest of this 
reaction comes from the fact that the scale dependence of the cross section is well 
under control, and therefore, the theoretical predictions are reliable. We quote here 
one result of this paper, referring the interested reader to the original publication 
[tj. Fig. 15 displays the cross sections dcr/dxLL, where xll = p]_ (e -T?7 +e _T? 3 ' et )/2f5 7 , 
for various cuts on the rapidities 7/ 7 and rf et (E 1 is the energy of the initial photon). 
In the forward region where the rapidities r? 7 , r]^ et are large, the cross section is 
dominated by the resolved contribution. For some cuts, half of the cross section is 
due to the gluon distribution in the photon. But the observable range in x is small 
(x < .2) and the cross sections fairly small. However this type of data could be used 
to constrain the poorly known G 7 (x, Q 2 ). 
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Figure 15: The enhancement of the gluon contribution to the reaction 7 + p 
7 + jet + X due to cuts on the rapidities ry 7 and r/j et . 













































7 Conclusions 


We have proposed a new set of next-to-leading order parton distributions for real 
photons. We work in the A4S scheme with a variable number of massless flavors, 
keeping however the massive correction terms for charm and beauty quarks, in the 
direct contribution to the J 7 ] structure function. The perturbative contribution is 
assumed to vanish below a Ql value and a VDM type parametrisation is used for the 
non-perturbative input at this value. We give a detailed discussion on how to isolate 
the scheme invariant piece of the structure function to which the physical VDM input 
should be applied. The distributions are set up such that various input parameters 
can be easily changed: these are the value of Qq, the overall normalisation of the non- 
perturbative input and the x-dependence of the gluon which is poorly constrained by 
the photon structure function data. Using the LEP data, an attempt to decorrelate 
the perturbative contribution, controlled by Q q, from the non-perturbative one is 
made. It turns out that the error bars are too large to determine unambiguously 
the value of Q q. A best fit to low and medium-Q 2 LEP data is then done which is 
also shown to yield a good agreement with the world data. A large dispersion in 
the value of the best fit parameters is observed when a minimum y 2 fit is performed 
for each set of LEP data independently. We therefore propose several sets of parton 
distributions to “frame” the data. Since the gluon distribution cannot be constrained 
from the deep inelastic data, we suggest to consider the photoproduction of photons, 
hadrons or jets at large transverse momentum. 
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A - Appendix A 


In this appendix, we give a derivation of expression (EH). We start from the 
reaction in which the target photon, instead of being real, has a small virtuality 
p 2 (— p 2 = P 2 -C Q 2 ), but large enough for the perturbative approach to be valid. 
This allows us to study the structure of the HO corrections to P] and to under¬ 
stand how to take the real photon limit p 2 —> 0. The transverse structure function 
JF^(Q 2 , P 2 , u) (transverse with respect to the polarisation of the target photon) can 
be written 


(Q\P\n) = C» (A.1) 

1 

where, for simplicity, we only consider the Non Singlet contribution. 

This expression is obtained in resumming all the ln(Q 2 /P 2 ) in the quark distribu¬ 
tion qf s (Q 2 , -P 2 , n), whereas C q (n ) and C^ s (n) are expansions in power of a s (Q 2 ). 
This procedure defines a factorization scheme called “Virtual Factorization Scheme” 
in ref. Pi- In this scheme the Wilson coefficient C q (n) and the direct term C^ s (n) 
are known since the work of Uematsu and Walsh (2] • 

The quark distribution is a solution of eq. (Q (we drop the index NS and the 
moment variable n ) 


o o ra s (Q 2 ) d\ p m 

Vf(Q ,p ) — a f , Tim KW e x ^ qq{ ) 

J otc' T ~‘ ON 


(A.2) 


s(P 2 ) /3(A) 

To simplify the notation further, we consider one quark species and drop the 
charge factors which are present in <jfk q (\) of (1A.2I) (that we shall note k( A)) and 
in C^ s . With this convention, the direct term is written 


C^(x) = — 6{ \(x 2 + (1 — x) 2 ) In—|- 2x(l — x) — 1 
27t L x 

+ ( x 2 + (1 — x ) 2 ^ In — + 6x(l — x) — l| . 


(A.3) 
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It is worth noting the following points. First C 7 (x) is scheme dependent and 
different from the MS expression (EZD- However it is easy to move from the Virtual 
Scheme to the MS scheme by keeping in mind that expression (EIH must be scheme 
invariant, which leads to 


a_ f kW MS (n) _ fcW(n) \ 
2^r [ pfi\ n ) P${n)) 
where, from (12.71) and (A.3) 


Cf 5 (n) - C» = AC'(n) 


(A.4) 


A C(x) = ^-6 {( x 2 + (1 — a;) 2 ) In x(l — x) + lj (A.5) 

in the x-space. 

Second, the direct terms C 1 are target dependent and depend on the regulariza¬ 
tion used to avoid a collinear divergence in the calculation of the box diagram. For 
instance, in dimensional regularization, the target dependent part has been calcu¬ 
lated in section 3, and is given by Cl /c (13.41) (in which ej is dropped). For a virtual 
photon, we obtain the first line of COJ) hu, whereas the second line is universal 
and equal to the equivalent MS expression. Therefore we can write in general 


C^x) = C?(x)+C”{x) . (A. 6 ) 

With the notation “col”, we indicate that the target dependent part of C 7 comes 
from the lower limit of the integration over k 2 (cf. expression 13.21 and 13 . 31) . A similar 
decomposition exists for Ad 1 ) 


Jfed) = kW** + Ad 1)C7 (A.7) 

because the target dependent terms, which are present in C 7 (x), also appear in the 
course of the calculation of A;*A under the form C^ ol {n) P^\n) (this point has been 
discussed in detail in ref. (521). Therefore the combination 
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“ <:, 1 ,(n) + c» = -^^ + c» 


2tc p/°) 


(A.8) 


qq yU) 
2 r>2 


2 vr p£\n) 

and consequently JF 2 (Q 2 , P 2 , n) (when Q 2 is very large) does not depend on the 
details of the target. 

We are now ready to study the limit P 2 — > 0 of expression (IA.2j) . With this 
aim in view, we introduce an intermediate scale Q q which allows us to isolate the 
P 2 -dependent part of El 


1(Q 2 ,P 2 ) = <,(Q 2 ,QD + £r°!mt) 


J. 


“s(Q ) dX f 

WP 


PqqP') 


, ( p.) m - ’ lA,9) 

where q(Q 2 , Qq) is given by (IA.2I) in which P 2 is replaced by Ql (in section 3, Q q is 
the limit between the perturbative and the non-perturbative domains). Using the 
notation 


fc(A) = — (fcW + — k^ U + A k n~ l ] = k u (X) + — — 

V ; 2vr \ 2t r 2vr j V ; 2tt 2tt 


(A.10) 


we rewrite the integral in (A.9) 


fQo dk 2 , n ... r v 5 dPlp a 

/ -—k u (X) IprPqi - 

JP 2 fc 2 27T 


Qq dk ' 2 ■ 


a 

pi Y° l 

t 

1 - 
l 

(«.(Q§r 

p(o) \ 

~ 2 ~W 


2tt 

p(0) 

JTqq 

wn, 

) 

> 


r a (Q ) dx 1 p (\/\ 
( .gh(g2) w7) ^ j 


f A.l 1) 

After extracting from (IA.11I) the P 2 -independent, but target dependent term 
k {1)C ° l /P£\ we define 


H(Ql,P 2 ) = 

and rewrite Pf as 


[ QU 4k«el 

Ip 2 k 2 


Q o dk 12 p a 

k 2 UT^qq -|- 


^ (i) (<*.m 

2p p<;> U,(p), 


p(0) 

„. , 2 qq 

\ 2 Po 


(A.12) 
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T. 


\Q\ P\ n) = C, q(Q\ Q\) + C, [H(Q„, P 2 ) 


qcoI 


S' 

o 


!s(Q 2 ) dX p 
s(Q 2 ) M 


+ cy + c: 


col 


(A-13) 

/■-i \coZ 

because = ^r—(or- Let us now consider the limit P 2 —> 0 with Ql being the 

scale below which the perturbative approach has no meaning. As a consequence, 

the perturbative expression (IA.12D of H{Q q, P 2 ) is no longer valid and we can only 

say that H(Qq, P 2 ) contains all the non-perturbative contributions needed to define 

P 2 in the real limit. Now we recognize in (I A. 1 Mil the structure of the input proposed 

in formula (|3. 71) if we identify H(Ql) = lim H(Ql,P 2 ) to q NP (Ql) + q NP (Ql). 

p 2 —> o 

Expression (1A.13I) has been established in the virtual factorization scheme. But 
we can easily obtain a similar expression in the MS scheme by writing fcL) — kW MS -f 
Afc^ 1 ) in the expression (EH for q(Q 2 , Qq). This change generates a MS distribution 
q MS (Q 2 ,Qo) and a term (cf. (12.121)1 


A/c (1 ) 

p(0) 

Pqq 


'otsiQ 2 


-2P$/Po' 


®s{Ql 


(A.14) 


which is combined with CZ° l 1 


_ ( a s (Q 2 ) \ 

WQijL/ 


2^^-2P 9 ( 9 0 ) //3o N 


of (I A. 131) to give 


c: 


col,MS 
7 


1 - 


V 


«s(Q 2 

K a s (Qo 


-2P$/0o' 


so that (IA. 131) can be written in terms of MS perturbative expressions, together with 
a non-perturbative contribution H(Q q) which remains invariant under this change. 
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B - Appendix B 


The parametrizations discussed in this paper are available in the form of a FOR¬ 
TRAN code allowing the users to select the parton distributions they are interested 
in. The conventions we use are described in table 1 which displays abbreviations of 
the general notation AFG04 (Qq, C np ,Pio)- 



Ql = .34 GeV 2 

Ql = .70 GeV 2 

Ql = .97 GeV 2 

Pio = 1-0 
(hard gluon) 


AFG04_BF_1.0 
(C np = .78) 


Pio = 1-9 
(default) 

AFG04_LW 
{C np = .60) 

AFG04_BF 
(C np = .78) 

AFG04_HG 
(C np = .84) 

Pio = 4.0 
(soft gluon) 


AFG04_BF_4.0 
(C np = .78) 



Table 1 


The values of C np indicated in brackets are the default values ; the users can choose 
other values (after carefully reading section 5). These parametrizations can be down¬ 
loaded from the site http://www.lapp.in2p3.fr/lapth/PHOX_FAMILY/main.htmh 
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